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Abstract. Recently, Liang & Zhang (2005) found a tight correlation involving only observable quantities, namely the isotropic 
emitted energy £7 7l i 80 , the energy of the peak of the prompt spectrum Ep, and the jet break time if*. This phenomenological 
correlation can have a first explanation in the framework of jetted fireballs, whose semiaperture angle 9j is indeed measured 
by the jet break time t' y By correcting -E 7 ,i so for the angle 9j one obtains the so called Ghirlanda correlation which links the 
collimation corrected energy E-y and E' v , There are two ways to derive 8j from t\ in the "standard" scenario, corresponding 
to an homogeneous or instead to a wind-like circumburst medium. We show that the Ghirlanda correlation with a wind-like 
circumburst medium is as tight as (if not tighter) than the Ghirlanda correlation found in the case of an homogeneous medium. 
There are therefore two Ghirlanda correlations, both entirely consistent with the phenomenological Liang & Zhang relation. We 
then suggest to consider the difference between the observed correlations and the ones one would see in the comoving frame 
(i.e. moving with the same bulk Lorentz factor of the fireball). Since both E p and E 1 transform in the same way, the wind-like 
Ghirlanda relation, which is linear, remains linear also in the comoving frame, no matter the distribution of bulk Lorentz factors. 
Instead, in the homogeneous density case, one is forced to assume the existence of a strict relation between the bulk Lorentz 
factor and the total energy, which in turn put constraints on the radiation mechanisms of the prompt emission. The wind-like 
Ghirlanda correlation, being linear, corresponds to different bursts having the same number of photons. 
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1. Introduction 

The correlation between the collimation corrected emitted en- 
ergy in the prompt phase of GRBs (-E 7 ) and the peak energy 
of the vE v prompt spectrum (E p ), the so-called "Ghirlanda" 
relation, is the fundamental tool for the cosmological use of 
GRBs. 

Its importance calls for a robust and convincing expla- 
nation. Possible ideas and interpretations have already been 
proposed, but requiring some ad hoc assumptions: differ- 
ent authors (Eichler & Levinson 2004; Levinson & Eichler 
2005; Yamazaki, Ioka & Nakamura 2004; Toma, Yamazaki & 
Nakamura 2005) have underlined the importance of viewing 
angle effects assuming different geometries of the fireball (an- 
nular or patchy); while Rees & Meszaros (2005) pointed out 
that a strict relation between total energy and typical peak fre- 
quency can be understood in a easier way if the underlying 
emission process is thermal, therefore suggesting that dissipa- 
tive processes in the photosphere of the fireball can increase the 
thermal (black-body like) photon content of the fireball itself. 
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The collimation corrected energy E 1 is derived by multi- 
plying the isotropic emitted energy -E 7 .i SO by the factor (1 — 
cos#j), where 6*j is the semiaperture angle of the jet. To de- 
rive it, one must assume the uniform jet model and the density 
profile of the circumburst medium (e.g. homogeneous or wind- 
like). 

The Ghirlanda correlation has been derived for an homoge- 
neous density (i.e. n constant) and there are a few critical points 
concerning its derivation which are now under discussion. 

Firstly, in this scenario the jet opening angle is (Sari et al. 
1999) 

0J = 0.161 (J^) 3/8 (-^-) 1/S (1) 

where z is the redshift, r; 7 is the radiative efficiency and t- h a 
is the achromatic break time, measured in days, of the after- 
glow lightcurve . The efficiency r] 1 relates the isotropic ki- 
netic energy of the fireball after the prompt phase, Sk,iso> to 
the prompt emitted energy E lt i SO , through E'k.iso = E 1: i SO /rj 1 . 
This implicitly assumes that rj 1 -C 1 otherwise the remaining 

1 Here we adopt the notation Q — 10 X Q X , and use cgs units unless 
otherwise noted. 
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kinetic energy after the prompt emission is instead -Ek.iso = 
^7,iso(l — Vj)/ r li- This efficiency, in principle, could be dif- 
ferent from bursts to bursts, but in the absence of any hints of 
how its value changes as a function of other properties of the 
bursts and favoured by its low power in Eq.l, one assumes a 
constant value for all bursts, i.e. ry 7 = 0.2 (after its first use 
by Frail et al. 2001, following the estimate of this parameter in 
GRB 970508). 

Secondly, the density of the circumburst medium n can in 
principle be estimated through accurate fits to the lightcurves 
at different frequencies [or, equivalently, fitting the spectral en- 
ergy distributions (SEDs) at different times], but the fact that 
the emitted synchrotron spectrum is insensitive to n in the 
regime of fast cooling makes the n estimates somewhat un- 
certain. Furthermore, only for a minority of bursts we have 
enough data to constrain n even in a rather poor way (see e.g. 
Panaitescu & Kumar 2000, 2001). From this partial informa- 
tion, however, the estimated values of n range from ~ 0.1 to 
~ 30 cm -3 with a clear preference for an homogeneous den- 
sity scenario. Wind density profiles are acceptable and even 
preferred in few cases, in which however one cannot exclude 
the homogeneous density case. This led to the choice, made 
by Ghirlanda, Ghisellini & Lazzati (2004, hereafter GGL04) 
to assume an homogeneous density scenario for all bursts and 
assign, to all cases in which n was not estimated, a range of 
possible values of n, from 1 to 10 cm~ 3 , to derive the value 
of 6j using Eq.^ and, more importantly, the associated error. 
However, if GRBs are expected to originate from the death of 
very massive stars, the wind density profile appears as the most 
natural outcome of the final stages of the evolution of the burst 
progenitor. 

These issues on one hand must caution us about the pos- 
sible scatter of the points in the Ghirlanda correlation, which 
may well be the results of bursts having a distribution (instead 
of a single value) of n and r\ 1 , but on the other hand the ex- 
tremely small scatter found suggests that, for the considered 
bursts, these values are indeed clustered in a small range. 

In any case, these concerns (i.e. a possible distribution of n 
and ?/ 7 values) have been completely overcome by the finding, 
by Liang & Zhang (2005, hereafter LZ05), of a phenomeno- 
logical and model-independent correlation between E lt i so , E' p 
and t'y By considering 15 GRBs, and a flat Universe cosmol- 
ogy with £1m = 0.28 and ho — 0.713, the correlation found by 
LZ05 takes the form: 



£ 7 ,iso, 52 = (0.85 ±0.21) 



EL 



100 keV 



1.94±0.17 



/-1.24±0.23 



j-d 



(2) 

where primed quantities are calculated in the rest frame of the 
GRB, i.e. H = t/(l + z) and E' p = E p (l + z). The scatter of 
the data points around this correlation is small enough to en- 
able LZ05 to use it for finding constraints to the cosmological 
parameters. 

The main aim of the present paper is to discuss a few steps 
which we think necessary in order to deepen our understanding 
of the spectral/energy correlations in GRB, even if we do not 
claim to arrive to a complete or satisfactory interpretation. The 
Ghirlanda correlation (although derived earlier than the LZ05 



relation) should be considered as a first step towards the expla- 
nation of the purely phenomenological LZ05 relation. 

The first step is to demonstrate that the LZ05 correlation is 
equivalent to the Ghirlanda correlation. Secondly, we explore 
what happens to the Ghirlanda correlation if, instead of an ho- 
mogeneous medium, we assume that the density is distributed 
with a r~ 2 wind profile. We demonstrate that the LZ05 is again 
consistent with this new Ghirlanda-wind correlation. 

We then have not one, but two possibilities of relating E'p 
and E 1 , both consistent with the LZ05 correlation. We argue 
that the new (wind-like) Ghirlanda correlation cannot be easily 
discarded on the basis of the fit to the afterglow SEDs, which 
prefer the uniform density case. The present small sample of 
GRBs does not allow to test the two possibilities, but we men- 
tion one test to do so when estimates of the initial bulk Lorentz 
factor will be available. 

We then discuss a third step, pointing out that the observed 
E p and E 1 are affected by the relativistic motion of the fire- 
ball. In the simplest and standard scenario, which assumes that 
the observer's line of sight is within the jet opening angle and 
that the jet is homogeneous, E p and E~ are both boosted by a 
factor ~ 2r, where T is the bulk Lorentz factor of the fireball. 
Therefore what we see is an apparent correlation. Imposing 
that some correlation (even different from the observed one) 
exists in the comoving frame allows to put interesting con- 
straints on both the dynamics and the emission processes of 
bursts. 

As part of our effort, we report the updated tables for the 
relevant parameters used in deriving all the results and correla- 
tions discussed in the present paper. We discuss case-by-case 
the differences and changes of these quantities with respect to 
previously published papers. 

2. The sample 

LZ05 found their correlation (Eq.|2j by using 15 bursts. Most 
of them are the same bursts used by GGL04, but there are 
some differences in the reported parameters for the common 
GRBs. Also, since the publication of GGL04 some (primarily 
spectral) parameters have been updated (as a consequence of 
refined analysis) and published in the literature. Finally, other 
bursts have been detected recently, which bring the total num- 
ber of GRBs with "useful" data (i.e. z, E p and tj, excluding up- 
per/lower limits) to 1 8 (at the time of writing, September 2005). 
We stress that our sample includes only those bursts with se- 
cure measurements of z, E pca ^, tj. Consistency checks of the 
upper/lower limits can be performed (as shown in GGL04). For 
these reasons we present in Tab^our sample of 18 GRBs with 
the relevant input data and corresponding references which rep- 
resents the most updated collection of published parameters 
which are relevant for our analysis. We detail in the follow- 
ing all the relevant differences of this table with respect to that 
reported in GGL04 and LZ05. 

Tab.^contains all GRBs reported by LZ05 with the excep- 
tion of GRB 021211, and with the addition of 4 GRBs (GRB 
970828, GRB 990705, GRB 041006, GRB 050525). Our ex- 
clusion of GRB 021211 is motivated by the extremely uncer- 
tain determination of tj due to the scarcity of afterglow data 
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Table 1. Input parameters for the bursts of our sample, a and (3 are the photon spectral indices of the prompt emission spectrum 
and E p represents the (observed) peak energy of the vF v spectrum. When the errors are not given in the original reference, 
we assumed an average error (values in square brackets), otherwise we list the originally reported error (in round brackets). 
"References are given in order for the redshift (z) and for the spectral parameters (a, (3, fluence and its energy interval): 1 
Djorgovsky et al. 2001; 2 Jimenez et al. 2001; 3 Djorgovsky et al. 1998; 4 Hjorth et al. 1999; 5 Amati et al. 2002; 6 Vreeswijk et 
al. 2001; 7 Amati et al. 2000; 8 Vreeswijk et al. 1999; 9 Amati 2004; 10 Hjorth et al. 2003; 1 1 Atteia et al. 2005; 12 Price et al. 
2003; 13 Barth et al. 2003; 14 Barraud et al. 2003; 15 Moller et al. 2002; 16 Sakamoto et al. 2005; 17 Greiner et al. 200 3a; 18 Rol 
et al. 2003; 19 Greiner et al. 2003b; 20 Vanderspek et al. 2004; 21 Weidinger et al. 2003; 22 http://space.mit.edu/HETE/Bursts/; 
23 Foley et al. 2005; 24 Blustin et al. 2005; ^References are given in order for the observed jet break time ij and for the density 
n when present: 25 Bloom et al. 2003; 26 Kulkarni et al. 1999; 27 Israel et al. 1999; 28 Bjornsson et al. 2001; 29 Jakobsson et 
al. 2003; 30 Berger et al. 2002; 31 Holland et al. 2003; 32 Klose et al. 2004; 33 Andersen et al. 2003; 34 Berger et al. 2003; 35 
Jakobsson et al. 2004; 36 Panaitescu & Kumar 2002; 37 Stanek et al. 2005; 38 Frail et al. 2005. 39 Blustin et al. 2005 



and to the likely "contamination" of a supernova (Delia Valle 
et al. 2003). The value of ij reported by LZ05 (1.4 day, see also 
Holland et al. 2004) should be considered only as a lower limit 
to tj (as assumed in GGL04). 

For GRBs listed in Tab. ^ which are also present in LZ04 
we have the following (minor) differences: 

- GRB 011211: the value of tj reported by Jakobsson et al. 
(2003) is <j = 1.56 ± 0.02 days. This has been derived 
by fitting the light curve with a broken power law, which 
very likely severely underestimates the error on this quan- 
tity. Note that the \ 2 of the original fit is large. Differently 
from LZ05, who used the above error, we have set the er- 
ror to 10% (i.e. the average error on the observed jet break 
time) of the value of t- 3 . 

- GRB 020124: the spectral parameters have been updated 
by Atteia et al. (2005) with respect to the ones assumed by 
LZ05 from the Sakamoto et al. (2005). 

- GRB 020405: for the value of E p listed in Tab.Q]we assume 
that in the original reference (Price et al. 2003), the authors 



presented the value Eq of the Band spectrum. This corre- 
sponds to the e-folding energy of the exponential rollover 
of this model and it allows to derive the peak energy as 
E p = Eg(a + 2) (where a is the photon spectral index of 
the low energy power law component of the Band model). 
Note also that if [3 — —1.87 ± 0.23, as reported in Price et 
al. 2003, the spectral peak energy is unconstrained within 
the observational energy band. However, the reported un- 
certainty (0.23) makes the high energy spectral component 
consistent with (3 < —2. For this reason we assumed the 
lowest value, i.e./3 = —2.1, to compute the isotropic equiv- 
alent energy of this burst. 
- GRB 020813: we take the spectral parameters of this GRB 
from Barraud et al. (2003), while LZ05 use the values re- 
ported in the Sakamoto et al. (2005). The reason of our 
choice is in the fact that, with the [3 = —1.57 ± 0.04 re- 
ported by Sakamoto et al. (2005), one could not define the 
value of E p = 140 ± 14 keV (which instead requires at 
least (3 < —2), which is reported in Sakamoto et al. (2005). 
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- GRB 021004: the value of tj = 4.74 days comes from 
Holland et al. (2003), and it is the same (as well as the 
quoted reference) of LZ05. What is different is the error: 
Holland et al. (2003) report t } = 4.74+°;^ days. LZ05 take 
0.14 days as the error in this quantity, while we take 0.5 
days to better approximate the (asymmetric) error. 

- GRB 030329: We have updated the spectral parameters of 
this bursts according to the published paper by Vanderspek 
et al. 2004 (which are somewhat different from those pre- 
sented in the preprint version of the same article). 

- GRB 050525: the prompt emission and the early afterglow 
detected by the Swift BAT, UVOT and XRT instruments 
have been recently analyzed by Blustin et al. (2005). In this 
paper the authors suggest the presence of a jet break in the 
very early afterglow lightcurve, i.e. tj =0.16 or 0.2 days 
according to two different model fit (the second value is 
found including also the data of Klotz et al. 2005, while the 
first value refers to the Swift data only). On the other hand, 
Mirabal et al. (2005) noted a jet break time at 0.4 days af- 
ter trigger, based on a large collection of data (120 frames, 
still unpublished) taken with the 2.4 MDM meter telescope. 
Note also that the fit leading to tj = 0.2 days assumes that 
the optical afterglow has a "jump" in its flux, but that there- 
after continues to decay normally. In other words, this "dis- 
continuity" is not treated as a "bump" in the lightcurve (as 
often seen in other bursts) which would have implied that 
the afterglow is composed by two contributions (the normal 
afterglow plus the re-brightening component). Fitting with 
this model would result in a larger value of tj . While wait- 
ing for a joint fit with all the data available, and also with 
different models, we decided to use tj ~ 0.28 ± 0.12 days: 
this value is intermediate between the Blustin et al. values 
and the Mirabal et al. value, with an error which encom- 
passes all values. 

In Tab.|2]we report the values of the rest frame peak energy 
E' and of the isotropic equivalent energy E ly - lso . The values of 
-E 7j ; so have been taken directly from Amati et al. (2002) in the 
case of GRBs detected by BeppoSax and listed in that paper, 
but we converted these values to our cosmology (i.e. we use 
h = 0.7 while Amati et al. 2002 used h = 0.65). 

Tab.Q an d Tab.|2]are the updated version of the tables pre- 
sented in GGL04 which were composed by 15 "usable" GRBs. 
With respect to that paper, our present sample of 18 GRBs 
comprises 3 new GRBs (GRB 021004, GRB 041006 and GRB 
050525) for which z, E p and tj have been published and some 
of the parameters have changed due to updates which appeared 
in the literature. 

- GRB 991216: the value of the density n = 4.70±^ j was 
estimated by Panaitescu & Kumar (2002). We have here 
assumed a symmetric error equal to the logarithmic aver- 
age, and taken n = 4.7 ± 3.5. This differs slightly from 
what assumed in GGL04 (n = 4.7 ± 2.3). 

- GRB 1 1 2 1 1 : we changed the reference for the spectral pa- 
rameters, which is Amati (2004), and not Amati et al. 2002 
as given in GGL04. We recalculated the value of i? 7 .iso 
with h = 0.7. 



- GRB 020124: we updated the spectral parameters, now 
taken from Atteia et al. (2005). The main difference con- 
cerns E p = 120 keV, instead of the value of 93 keV re- 
ported in GGL04. 

- GRB 030226 and GRB 030328: we now use the spectral 
parameters reported in Sakamoto et al. (2005) (instead of 
Atteia 2003) and the recalculated E ly i so . 

- GRB 030329: as mentioned above, we updated the spectral 
parameters according to the published paper by Vanderspek 
et al. 2004. The listed value of n — 2.2 cm~ 3 comes from 
Frail et al. 2005 (GGL04, instead, quoted n = 1 cm -3 ): 
the associated error encompasses the three possible values 
listed in Tab. 2 of Frail et al. 2005. 

- GRB 030429: we updated the spectral parameters with 
Sakamoto et al. (2005) instead of those assumed in GGL04 
which were taken from the Hete-2 web page. 



3. The spectral-energy correlations 

With the updated sample of 18 GRBs reported in Tab.^we first 
refit the empirical LZ05 correlation among E ly - lso , E' p and t' y 
We also give the updated version of the Ghirlanda correlation 
(GGL04 and Ghirlanda et al. 2005) in the case of a homoge- 
neous density profile. Finally we present the Ghirlanda corre- 
lation in the case of a wind density profile. 

3. 1. The Liang-Zhang correlation revisited 

The method used by LZ05 to find the correlation between 
E lt i so , E' p and tj is a multivariate linear regression. The signifi- 
cance of the multivariate regression is estimated through the F- 
test and through the Spearman r s coefficient between log E ly i so 
calculated through Eq.[2]and the same quantity directly calcu- 
lated through 



E 



7, iso 



4ird 2 L S 7 k 
l + z 



(3) 



where is the luminosity distance, 5 7 the 7-ray fluence in 
the observed energy band, k is the bolometric correction factor 
needed to find the energy emitted in a fixed energy range (here, 
1-10000 keV) in the rest frame of the source. 

We have used a different method, which enable us to weight 
the multidimensional fit for the errors on the three independent 
variables Ej^ so , E' p and tj. By extending to the three dimen- 
sional space the procedure for the fit of a straight line to data 
with errors on two coordinates (Press et al. 1999) we use the 
X 2 statistics to find the best fit. 

Assuming a 57m = 0.3 and h = £1a = 0.7 cosmology, we 

find 



E. 



7, iso, 53 



(1.12±0.11) 



EL 



295 keV 



1.93±0.17 



t 



0.51d 



-1.08±0.17 



(4) 

with a reduced \ 2 = 149. The 3D plot of the data points in the 
i? 7 .iso, E' p and tj space and the best fit plane as defined by Eq.|4] 
are represented in Fig.^ Similarly to what has been done in 2D 
(GGL04) we can define the scatter of the data points around the 
best fit plane through their distance computed perpendicular to 
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GRB 


K 




9i 




i5 7 ,n=3 








keV 


erg 


deg 


erg 


erg 


deg 


erg 


970828 


583 [117] 


2.96e53 [0.35] 


5.91 ±0.79 


1.57e51 (0.46) 




3.40(0.19) 


5.21e50 (0.84) 


980703 


499 [100] 


6.9e52 (0.82) 


11.02 ± 0.80 


1.27e51 (0.24) 


7.29e50 (1.20) 


5.45 (0.26) 


3.12e50 (0.48) 


990123 


2031 (161) 


2.39e54 (0.28) 


3.98 ± 0.57 


5.76e51 (1.78) 




1.84 (0.12) 


1.24e51 (0.22) 


990510 


422 (42) 


1.78e53 [0.19] 


3.74 ± 0.28 


3.80e50 (0.69) 


6.80e50 (0.98) 


3.31 (0.14) 


2.98e50 (0.40) 


990705 


348 (28) 


1.82e53 (0.23) 


4.78 ± 0.66 


6.33e50 (1.92) 




3.20(0.19) 


2.84e50 (0.50) 


990712 


93 (16) 


6.72e51 (1.29) 


9.47 ± 1.20 


9.16e49 (2.90) 




8.75 (0.51) 


7.82e49 (1.76) 


991216 


642 [129] 


6.75e53 [0.81] 


4.44 ± 0.70 


2.03e51 (6.79) 


1.81e51 (0.50) 


2.36 (0.21) 


5.72e50 (1.24) 


011211 


185 (25) 


5.4e52 (0.6) 


5.38 ± 0.66 


2.38e50 (0.64) 




4.24(0.16) 


1.48e50 (0.20) 


020124 


390 [113] 


2.61e53 (0.18) 


5.07 ± 0.64 


1.02e51 (0.27) 




3.13(0.12) 


3.90e50 (0.40) 


020405 


617 [124] 


1.25e53 [0.13] 


6.27 ± 1.03 


7.48e50 (2.58) 




4.08 (0.34) 


3.17e50 (0.63) 


020813 


478 [95] 


5.78e53 [0.58] 


2.8 ± 0.36 


6.89e50 (1.88) 




1.85 (0.08) 


3.00e50 (0.40) 


021004 


267 (117) 


3.38e52 (0.78) 


8.47 ± 1.06 


3.69e50 (1.25) 




6.20 (0.40) 


1.98e50 (0.52) 


030226 


290 [63] 


5.43e52 (0.68) 


4.71 ±0.58 


1.84e50 (0.51) 




3.88 (0.17) 


1.24e50 (0.19) 


030328 


318 [33] 


3.68e53 [0.37] 


3.58 ± 0.45 


7.18e50 (1.93) 




2.35 (0.10) 


3.09e50 (0.40) 


030329 


82 (2) 


1.62e52 [0.16] 


5.69 ± 0.50 


7.99e49 (1.62) 


8.63e49 (1.56) 


5.52 (0.31) 


7.5 le49 (1.13) 


030429 


128 [26] 


1.42e52 (0.23) 


6.3 ± 1.52 


8.57e49 (4.37) 




5.88 (0.88) 


7.48e49 (2.55) 


041006 


108 [22] 


6.92e52 [0.7] 


2.79 ± 0.41 


8.18e49 (2.57) 




2.62(0.18) 


7.25e49 (1.23) 


050525 


135 [8] 


2.3e52 [0.3] 


4.04 ± 0.8 


5.73e49 (2.34) 




4.04 (0.45) 


5.72e49 (1.40) 



Table 2. The rest frame peak energy E' p , the collimation corrected energy -E 7 .i SO are calculated with the parameters reported in 
Tab.^ The semiaperture angle 9j is calculated with Eq.[^in the case of an homogeneous medium and the collimation corrected 
energy E 7 is reported. We also report the values of i? 7 . n= 3 assuming n = 3 cm~ 3 for the 4 GRBs which have a different estimate 
of n as reported in Tab.^ In the case of a wind medium the semiaperture angle 8j. w (calculated by Eq.[6} and the corresponding 
E~y. w are reported. When the errors in E p (i? 7i ; so ) are not given in the original reference, we assume the average error of 20% 
(11%) (values in square brackets), otherwise we list the originally reported error (in round brackets). 



Comparing Eq.|4]with the original result of LZ05 we obtain 
a value of the <j exponent closer to unity, but still consistent 
with the value of LZ05 (which was also consistent with unity, 
due to the relatively larger uncertainty). As explained in the 
next section, a value equal to unity is crucial to make the LZ05 
correlation and the Ghirlanda correlation mutually consistent. 

We show in Fig. [3] the values of -E 7 .i so calculated through 
Eq. [3] as a function of -E 7) i so calculated through the best fit 
of the correlation E 1 \ so (E' p: ij) found using our data and our 
method, and assuming a Vlyi — 0.3 and h = fl\ = 0.7 cos- 
mology. As can be seen, there is a very good agreement. 

It is interesting, in view of the discussion of the following 
sections, also to fit the E lt i so , E' p and correlation by forcing 
the slope of the t'j to be -1 . Fixing this slope we find E 1 ,i so t'^ oc 
£i- 91±01 , with a reduced X l = 1-51. 

3.2. The updated Ghirlanda correlation 

Using the same data listed in Tab. ^ we calculate the up- 
dated version of the Ghirlanda correlation. We find a Spearman 
correlation coefficient r s =0.93 with a chance probability 
P=2.3xl0~ 8 . We report in Fig. 0] the updated correlation with 
the 18 GRBs reported in Tab.Q 




-0.5 0.0 0.5 

3D scatter 

Fig. 2. Histogram of the scatter of the data points computed 
perpendicular to the best fit plane in the 3D space of Fig. \l\ 
The solid line represents the gaussian fit which has a a = 0.18. 

this plane. The histogram of the scatter is reported in Fig|2]and, 
when fitted with a gaussian it has a a = 0.18. 
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Fig. 1. Three dimensional representation of the E~ l! i so (E' p , ij) correlation. Data points are from Tab|2| The plane is the best fit to 
the data points as represented by Eq.[4]which has a reduced Xi — 1-49. The red points lies above the best fit plane while the grey 
stars are below the plane. The height of each point with respect to the best fit plane is also shown. 




E iso [erg] 

Fig. 3. -E 7 ,iso calculated through the best fit multidimensional 
correlation E ly i so (E' p ,t'j) (Eq. as a function of the same 
quantity calculated through Eq.[3] 

We fitted this correlation with a powerlaw model account- 
ing for the errors on both variables, i.e. E' p and E 7 (using the 



routine fitexy of Press et al. 1999). For J7m = 0.3, £1a = h = 
0.7 we find 

/ jpi \ / \ 0.69±0.04 

— =M = (2.79 ±0.15) ( ^—] 

VIOOkeVy V ' V2.72 x 10 50 erg/ 

(5) 

with a reduced Xi = 1-4 for 16 degrees of freedom. The errors 
on its slope and normalization are calculated in the "barycen- 
ter" of E' p and E Jt where the slope and normalization errors 
are uncorrected (Press et al. 1999). We also note that the sim- 
plest linear regression fit (i.e. without accounting for errors on 
the variables) gives a slope of 0.6 (dotted line in Fig.|4}. 

Fig lUshows the correlation for the 18 GRBs and its best 
fit represented by Eq. |5] (solid line). This updated correlation 
has a slope which is consistent with the original value found in 
GGL04. We also computed the scatter of the data points around 
this correlation (insert of Fig0}. This scatter is defined as the 
distance in the log E' — log E 1 plane of each data point from 
the best fit correlation. We find that if fitted with a gaussian its 
standard deviation is a = 0.1, i.e. lower than the value origi- 
nally found by GGL04. 

For completeness, we also computed the Amati correlation 
with the 18 GRBs reported in Tab.Q By weighting for the er- 
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Fig. 4. The updated Ghirlanda correlation between the rest frame spectral peak energy E' and the collimation corrected energy 
Ery as found with the 18 GRBs reported in Tab.^ The solid line represents the best fit powerlaw model obtained accounting for 
the errors on both coordinates (Eq. |5j which has a reduced x? = 1-4 (16 dof) and a slope of 0.69±0.04. We also show the fit 
obtained with the simplest linear regression, i.e. without accounting for the errors on the coordinates (dotted line, slope equal to 
0.6). The circled point represents GRB 990510 which alone contributes to the 27% of the total \ 2 of the fitted model. The names 
of the 18 GRBs are also reported. The shaded areas represent the regions corresponding to the 1, 2 and 3a scatter around the best 
fit correlation. The insert reports the distribution (hatched histogram) of the scatter of the data points computed perpendicularly 
to the best correlation (solid line in main plot) and its gaussian fit (solid line in the insert) which has a a = 0.1. 



rors on E p and E 7t i SO , we find a relatively poor fit with a re- 
duced Xr = 5.22 and a best fit correlation S 7iiso oc ££-57±o.02 

3.3. The Ghirlanda correlation in the case of a wind 
density profile 

If the external medium is distributed with anr -2 density profile 
the semiaperture angle of the jet is related to the achromatic jet 
break through (Chevalier & Li 2000): 

9, w = 0.2016 V4 (-^-) V4 (6) 



where we assume n(r) = Ar~ 2 and A* is the value of A {A = 
Mw/(47ru w ) = 5 x 10 n j4* g cm -1 ) when setting the mass 
loss rate due to the wind M w = 10~ 5 Mq yr _1 and the wind 
velocity v w — 10 3 km s — 1 , according to the Wolf-Ray et wind 
physical conditions. 

In the wind case we use Eq.[6]to correct the isotropic en- 
ergies -E 7 .iso by the factor (1 — cos0j iW ). Given the few still 
uncertain estimates of the A* parameter, we assume the typical 
value (i.e. A* = 1) for all bursts neglecting for the moment the 
possible uncertainty on this parameter. 
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Fig. 5. The Ghirlanda correlation in the case of a wind profile of the external medium density as found with the 18 GRBs 
reported in Tab.^(col. 5). The values of E 1 are reported in Tab.[2] The solid line represents the best fit powerlaw model obtained 
accounting for the errors on both coordinates (Eq.0 which has a reduced Xr — 1-125 (16 dof) and a slope of 1.09±0.06. We also 
show the fit obtained with the simplest linear regression, i.e. without accounting for the errors on the coordinates (dotted line), 
which has a slope of 0.92. The circled point represents GRB 030326 which is giving the largest contribution (23%) to the best 
fit reduced \ 2 - The shaded regions represent the 1, 2 and 3er scatter around the best fit correlation. The names of the 18 GRBs 
are indicated. The insert reports the distribution (hatched histogram) of the scatter of the data points computed perpendicularly 
to the best correlation (solid line in the main plot) and its gaussian fit (solid line in the insert) which has a a = 0.08. 



In the wind case we find a correlation between E' p and E~ ( 
with a Spearman rank correlation coefficient r s = 0.92 (P = 
6.9 x 10~ 8 ). The fit with a powerlaw model gives 



100 keV 



= (3.0 ±0.16) 



E~, 



2.2 x 10 50 erg 



1.03±0.06 



(7) 



with a reduced Xi = 1-13 for 16 dof (see Fig.|5}. Note that the 
exponent of this new relation is entirely consistent with unity. 
The scatter of the points (insert of Fig. [5} around the best fit 
correlation is fitted by a gaussian with a = 0.08. 



Since we have no knowledge of the uncertainty associated 
with the ijA* parameter entering in Eq.[6] we estimate that, for 
the assumed typical value rjA* — 0.2, an error a v A, < 20% 
does not dominate the fit of the correlation (i.e. the reduced \? 
is not much smaller than 1). 

Since we have assumed that all the A* values are equal and 
have no errors, the resulting Xr of the wind case should be com- 
pared with the the case of homogeneous density assuming all 
the n values equal (we set n = 3 cm~ 3 ) with no error (see 
Tab. |3. This case is shown in Fig. In this case we obtain 
Xr — 1.4 to be compared with the x r = 1-125 of the wind 
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Fig. 6. Jet opening angles calculated for a homogeneous den- 
sity profile (6j) and for a wind density profile (6j <w ) for the 18 
GRBs reported in Tab. 2 (col. 4 and 7 respectively). 

case. We conclude that the wind case gives a somewhat better 
Xr and a somewhat tighter (smaller scatter) and a steeper corre- 
lation than the correlation found in the case of an homogeneous 
medium. Therefore, even if the wind density profile is not fa- 
vored by the afterglow model we cannot discard it on the basis 
of the Ghirlanda relation. 

The jet opening angles calculated in the case of a homoge- 
neous medium or in the case of a wind density profile, for the 
sample of 18 GRBs, is reported in Fig. [6] The angle calculated 
in the wind case (Eq. |6) is sistematically smaller than in the 
homogeneous medium case (Eq.Q. 

4. Consistency of the empirical correlation with 
the model dependent correlations 

In this section we will demonstrate that the Ghirlanda correla- 
tion either assuming an homogeneous circumburst density or 
a wind density profile and the LZ05 correlation are mutually 
consistent. This allows us to make some interesting considera- 
tions on the scatter of the Ghirlanda correlations. The fact that 
the LZ05 and the Ghirlanda correlation are mutually consistent 
have been already pointed out by LZ05 (see also Xu 2005) in 
the case of an homogeneous medium. 

For the simple analytical demonstration, we will consider a 
generic form of the LZ05 correlation, namely: 

(8) 



■p p/A.fB 



4.1. Homogeneous density 

Adopting the standard fireball scenario, assuming a uniform 
jet and an homogeneous circumburst density distribution, the 
relation between £j and 0- 3 is given by Eq. [J Inserting it into 
Eq.[8]one obtains 

>(mhr B/(3 ~ B) 

(9) 



E 



E 1 OC £^l/(3-B) 06 C B+l)/(3- 
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Fig. 7. The Ghirlanda correlation in the case of a homogeneous 
external medium assuming a universal value for the density, i.e 
n=3.0 without uncertainty, for comparison with the wind case. 
The solid line represent the best fit powerlaw model obtained 
accounting for the errors on both coordinates which has a re- 
duced x 2 r = 1-4 (16 dof) and a slope of 0.69±0.04. 



where we have used the small angle approximation (1 — 
cos0j) oc 6?. Eq. prelates five variables. 

We have shown that the Ghirlanda correlation is character- 
ized by a scatter a ~ 0.1 (Sec. 3.2 and Fig.|4). If this scatter is 
entirely due to a dispersion of the nr; 7 values (and very likely 
also to the errors of measurements on the observables) and not 
to 6*j, we derive the condition B = — 1 and Eq. [9] reduces to 
a relation between E 1 and E' p . This value of B is consistent 
with the value found from the fit of the LZ05 correlation (i.e. 
B = -1.08 ± 0.17) with the 18 GRBs of our sample. If the 
exponent B = — 1, then the slope g of the Ghirlanda correla- 
tion (i.e. E' p oc E!>) is related to the exponent A of the LZ05 
correlation through g = 4/ (3 A). 

On the other extreme, if the scatter of the Ghirlanda correla- 
tion is completely due to 9j (and to the errors of measurements 
on the observables) we may still derive a range of allowance 
for the parameter B. We computed the standard deviation a of 



the distribution of 9- 



6(fl+l)/(3-fl) 



as a function of B. This is 



represented in Fig. [8] by the solid line. If we compare a with 
the scatter of the Ghirlanda correlation (in the homogeneous 
density case - solid horizontal line in Fig.[8j we can find quite 
shallow constraints on B E (—1.3, —0.7). 

Clearly, the intermediate case corresponds to both nry 7 and 
Oj contributing to the scatter. If we consider the range of B as 
found by fitting the LZ05 correlation (Eq.0, i.e. B = -1.08 ± 
0.17 (shaded region in Fig. [8} then the term 0?(-B+i)/(3 b) can 



contribute at most for the 
Ghirlanda correlation. 



70% of the total scatter of the 



10 



L. Nava et al.: Spectral-energy correlations in long GRBs 




Fig. 8. Standard deviation of the distribution of 9 h {B) where 
h(B) is given by Eq.|9]and Eq.^|for the homogeneous (solid 
line and filled points) and wind density profile (long dashed line 
and open circles) as a function of the parameter B. The hori- 
zontal lines represent the scatter of the Ghirlanda correlation 
in the two scenarios (solid and long-dashed line respectively). 
The shaded region represents the (la) uncertainty on the pa- 
rameter B (=-1.08±0.17) found through the fit of the LZ05 
correlation with the 18 GRBs. 

4.2. Wind density profiles 

If tj is related to the semiaperture angle of the jet according to 
Eq.[6] the LZ05 correlation implies: 

Ej . iso ef w ~ e, oc E'/ni-B) 02 (B+ i)/(i-B) (m ^_ B/(1 _ B) 

(10) 

For the same considerations reported in the previous section, 
we derive B = — 1 if only the term nr/j contributes to the scat- 
ter of the Ghirlanda correlation in the wind density case (i.e. 
er ~ 0.08, Sec. 3.3 and Fig|5Jl. In this case the relation between 
the Ghirlanda correlation and the LZ05 correlation implies that 
E' p oc Ej /A . 

If we consider that the scatter of the Ghirlanda correlation 
(long dashed horizontal line in Fig. [8} is entirely due to the 
term g^- 8 ^^ 1 s ' we can derive even shallower (with re- 
spect to the homogeneous case) constraints on the parameter 
B G (—1.45,-0.65) (Fig. [H] - long dashed curve). Again in 
the intermediate case, i.e. both nry 7 and 0j contributing to the 
scatter of the Ghirlanda correlation, if we assume the range of 
possible values of B resulting from the fit of the LZ05 corre- 
lation than the term 9^ contributes at most for the 

^60% to the scatter observed in the Ghirlanda correlation. 

5. The Ghirlanda correlation in the comoving 
frame 

The values of E' p and E-y we derive are the ones seen in the 
GRB rest frame but not in the fireball comoving frame. 



In the simplest and standard scenario, which assumes that 
the observer's line of sight is within the jet opening angle and 
that the jet is homogeneous, E p and E 1 are both boosted by a 
factor ~ 2T, i.e. the blueshift of the photons. 

Then the comoving E™ m = E 1 /{2Y) and E™ m = 
E' p /(2T). Here T is the bulk Lorentz factor of the fireball emit- 
ting the prompt radiation. 

We then exploit the fact that we see a very tight correlation 
(in the rest, but not comoving, frame) to pose some limits on the 
physics (i.e. the dynamics and the radiation process) of bursts. 

Let us assume that the most general Ghirlanda correlation 
described by a generic power law, i.e. E' p oc E^, which can 
then represents the correlation for both the homogeneous and 
wind density profiles. In the comoving frame, each point must 
be corrected by the bulk Lorentz factor of that burst. If these 
bulk Lorentz factors are uncorrelated with E 1 or E' p , then in 
the comoving frame the correlation is destroyed. It would then 
seem very strange that two quantities that are not correlated in 
the comoving frame appear to be correlated in the rest frame. 
We then are obliged to assume that T is a function of E 1 (or, 
equivalently, of £^ om ). To this aim, let assume a simple power 
law relation: 

r OC {E c ° m ) X ' ; r OC (11) 

The two relations above are equivalent. In the comoving frame 
we have: 

E c p om oc (E™ m ) X9 ~ x+9 (12) 

Note that: 

- If x = (i.e. all bursts have the same Lorentz factor) then 
the slope in the comoving and in the rest frame is the same; 

- If g = 1 (i.e. wind case), the correlation is linear also in the 
comoving frame; 

- The exponent x can be thought as determined by a partic- 
ular dynamical model of the fireball. Once x is fixed, then 
the exponent (xg — x + g) appearing in Eq.[21should be 
explained by the radiative process. Viceversa, if we have 
reason to fix, through a specific radiation model, the expo- 
nent of Eq.^l then we have information on the dynamics 
of the fireball. 

We stress the fact that for the wind case (i.e. g — 1) we 
obtain a linear relation, whose slope is therefore "Lorentz in- 
variant" . This of course would greatly help to explain the ex- 
istence of the Ghirlanda correlation, since one of the main pa- 
rameters, the bulk Lorentz factors, does not enter to determine 
it. In other words, in the case of the wind, the explanation of 
the correlation should be found in the radiation process only, 
independently of the dynamics. 

6. Discussion 

We have shown that the model-independent correlation re- 
cently found by LZ05 between the isotropic emitted energy 
-E^iso, the rest frame peak energy E' p and the jet break time 
£j, calculated in the rest frame, is equivalent to the Ghirlanda 
correlation. 
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We have also shown that an even tighter correlation is 
found assuming that the circumburst density is distributed 
with an r~ 2 wind profile. Remarkably enough, this Ghirlanda- 
wind correlation is linear, and this slope is independent of the 
Lorentz correction needed to find out the correlation in the co- 
moving frame. 

We are aware that the wind-like distribution of the circum- 
burst medium is not favored by the existing afterglow mod- 
eling, but the advantage of having a linear Ghirlanda corre- 
lation is so great to justify a deeper analysis, to see if there 
can be some neglected effects which might be able to hide the 
presence of the wind. This is however out of the scope of the 
present paper, and we defer this issue to future studies. Here we 
only mention that one of the main assumption of the afterglow 
modeling might be particularly suspect, namely the hypothe- 
sis that the equipartition parameters e# and e e (i.e. the fraction 
of the dissipated energy converted in the magnetic field and in 
the electron energy, respectively) are kept fixed during the en- 
tire afterglow phase (while they have very different values from 
burst to burst). For instance, in the case of 63 oc T~ A , the syn- 
chrotron cooling is enhanced at later times with respect to the 
case of a constant cb, and this makes the light curve of a fire- 
ball moving in a wind circumburst environment to mimic the 
evolution of a fireball expanding in a uniform medium with a 
constant e# (the two cases becomes almost indistinguishable 
for A = 2). Consider also that all the pieces of evidence we 
have up to now point towards a massive stellar progenitor of 
GRBs, and it is difficult to understand why there is no sign of 
winds around such massive stars. The fact that e# changes (as 
long as its evolution law does not change), does produces a 
different (with respect to a not-evolving cb) decay law of the 
afterglow flux, but with no breaks. Therefore in this case eb 
does not enter in the estimate of the jet opening angle (as long 
as we are in the adiabatic regime). 

One interesting developement (Ghirlanda et al. 2005a) of 
having found a somewhat tighter and steeper Ghirlanda corre- 
lation in the wind density case (with respect to the homogenous 
case) is to use it to constrain the cosmological prameters simi- 
larly to what already done through the Ghirlanda correlation in 
the homogeneous case (Ghirlanda et al. 2004a; Firmani et al. 
2005). 

As pointed out by Ghirlanda, Ghisellini & Firmani 
(2005), the existence of the Ghirlanda correlation explains the 
Amati correlation between E y i so and E' p (Amati et al. 2002, 
Ghirlanda, Ghisellini & Lazzati 2004), and in particular ex- 
plains why the Amati correlation has a much larger scatter than 
the Ghirlanda correlation. In fact if GRBs are characterized by 
a distribution of semiaperture angles for each value of E' p , then 
one sees a variety of -E 7 4 SO -values for each value of E' p . If the 
(observed) distribution of aperture angles turns out to have a 
preferred value, where it peaks, then this naturally produces a 
correlation in the E' p -Ej^ so plane. Furthermore, it is conceiv- 
able that the bursts with spectroscopically measured redshifts 
are the brightest, hence with the smallest aperture angles (for a 
given redshift), hence lying at the large E Jt i so end of the real 
distribution. We have evidence that this is just what is happen- 
ing (Ghirlanda, Ghisellini & Firmani 2005). 



Therefore the Amati correlation can be easily explained by 
assuming the existence of i) the Ghirlanda correlation and ii) 
the existence of a (peaked) distribution of jet aperture angles. 
Thus what remains to be explained is the Ghirlanda correla- 
tion itself which is relating the intrinsic collimation corrected 
quantities. We have pointed out that there might be a differ- 
ence between the apparent and the comoving Ghirlanda corre- 
lation, and from the theoretical point of view it is the comov- 
ing one that should be explained. As Rees & Meszaros (2005) 
pointed out, a tight relation between the peak energy of the 
spectrum and the total emitted energy reminds of a thermal pro- 
cess, where the peak energy is a measure of a temperature. It is 
in this direction that we plan to investigate in the future. 

Finally, let us comment, regarding the linear, wind-like, 
Ghirlanda correlation, which implies that the number of "rel- 
evant" photons (the ones with energies close to E p ) is con- 
stant in all bursts, and approximately equal to N y = 10 57 , a 
number (coincidentally?) close to the number of baryons in 
one solar mass. In the "standard" scenario, in which the pri- 
mary energy is injected close to the putative newly born black 
hole in a high entropy form, the requirement of a fixed num- 
ber of photons N 1 translates in the requirement that the prod- 
uct of the injected energy E and the typical size of the injec- 
tion region R is constant. This can be seen in a simple way 
by noting that aT i R 3 /(kT) ~ iV 7 , where T is the tempera- 
ture of the initial blackbody. Since E ~ kTN 7 , one arrives 
to ER oc N^ 3 =const. In the framework of the standard sce- 
nario, in the absence of additional injection of energy, the num- 
ber of photons is conserved during the acceleration and the 
coasting phase. Then, when the fireball becomes transparent, 
this blackbody-like component has a fixed number of photons. 
But bursts with different baryon loading would become trans- 
parent at different times, meaning that different fractions of the 
energy initially contained in the radiation field have survived to 
the conversion to the bulk kinetic energy. This means different 
efficiency factors r\ 1 . And this is contrary to one of the main 
assumptions leading to the construction of the Ghirlanda corre- 
lation itself, which has been derived assuming the same ?y 7 for 
all bursts. 

If some extra dissipation of kinetic energy occurs after 
the acceleration phase, before the transparency radius, and the 
main radiation process is Comptonization, then we can recon- 
vert part of the kinetic energy into radiation energy leaving 
the number of photons unaltered. The available time to do it 
will be proportional to the transparency radius, so we have the 
positive feedback that sources which becomes transparent later 
(i.e. more baryon loaded, and hence with a radiation content 
which will be less energetic), have more possibilities to be re- 
energized by this putative dissipation (Rees & Meszaros 2005). 
This positive feedback can narrow the range of radiative effi- 
ciencies. In any case, we face the problem to explain why bursts 
with the same number of photons have different total energies. 

Apart from these theoretical considerations, there might be 
an "observational" way to see if the real E 7 -E' p is wind-like 
or constant density-like. Assume in fact that there is indeed a 
correlation between E p and E y and suppose that it will be pos- 
sible to estimate the bulk Lorentz factor of the fireball during 
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the emission of the prompt (that should be roughly equal to the 
bulk Lorentz factor of the very early afterglow, and controls 
the time at which there is the peak flux of the afterglow itself). 
Preliminary attempts to estimate it have already been done 
by e.g. Soderberg & Ramirez-Ruiz (2002) for GRB 990123. 
Then, if by "de-beaming" the E 1 and E' p values (i.e. comput- 
ing them in the comoving frame) in the case of an homoge- 
neous density one finds an equally tight correlation, this would 
strongly point in favor of the validity of the homogeneous den- 
sity hypothesis. Perhaps more importantly, one will also find 
the relation between V and E^° m . On the contrary, if the de- 
beamed quantities (calculated in the case of an homogeneous 
density) do not correlate, this will strongly argue in favor of the 
wind hypothesis (in such a case the de-beamed quantities still 
correlate). 
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